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Abstract The current generation of D-Wave quantum annealing processor is 
designed to minimize the energy of an Ising spin configuration whose pairwise 
interactions lie on the edges of a Chimera graph Cm,n,l- In order to solve 
an Ising spin problem with arbitrary pairwise interaction structure, the corre¬ 
sponding graph must be minor-embedded into a Chimera graph. We dehne a 
combinatorial class of native clique minors in Chimera graphs with vertex im¬ 
ages of uniform, near minimal size, and provide a polynomial-time algorithm 
that finds a maximum native clique minor in a given induced subgraph of a 
Chimera graph. These minors allow improvement over recent work and have 
immediate practical applications in the held of quantum annealing. 

Keywords Graph minor, clique minor. Chimera, graph embedding, adiabatic 
quantum computing, quantum annealing. 


1 Introduction and motivation 

D-Wave quantum annealing processors are designed to sample low-energy spin 
conhgurations in the Ising model using open-system quantum annealing m 
IH1I5] . Input to the processor consists of an Ising Hamiltonian (h,J), where 
h € K" is a vector of local fields and J G is a matrix of couplings, 

which we assume here to be symmetric. The energy of a spin conhguration 
s G {—1,1}" is given as 

E{s) = E{h,J,s)=s^Js + s^h. (1) 
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The output of an anneal (i.e. a run) of the processor is a low-energy state s, 
which consists of an Ising spin (either —1 or 1) for each qubit. 

Nonzero terms of J are physically realized using programmable couplers. 
These couplers only exist between certain pairs of physically proximate qubits. 
The input (/i, J) is therefore restricted such that if Jij ^ 0, there must be a 
coupler between qubit i and qubit j. We rephrase this in graph-theoretic terms: 
The connectivity graph of (h^J), which we denote by Gj, is the undirected 
graph on n vertices whose adjacency matrix has the same nonzero entries as 
J. Likewise, each processor has a hardware graph G representing the available 
qubits and couplers in the processor. For {h, J) to be input directly to the 
processor, Gj must be a subgraph of G. If this is not the case, we can in¬ 
directly input (h, J) to the hardware by embedding Gj as a graph minor of 
G. Implementing graph minors in the Ising model involves putting a strong 
ferromagnetic coupling Ji j 0 between any two adjacent vertices i,j of G 
in the same vertex image. This coupling compels multiple qubits to take the 
same spin, thus acting like a single logical qubit. The method is studied in 
greater detail elsewhere [dl lT^fTTl f^l^ . 

In this paper we consider the problem of finding large clique minors in 
the hardware graph. This is sufficient for minor-embedding any problem of 
appropriate size in a given hardware graph, and allows the study of random, 
fully-connected spin glass problems, as in recent work m- Klymko et al. first 
provided a polynomial-time algorithm for generating large clique minors in 
subgraphs of hardware graphs [10] . In Section we provide evidence that 
our algorithm uses fewer physical qubits and allows the embedding of larger 
minors. 


1.1 The Chimera graph and triangle embeddings 

D-Wave processors currently operate using a Chimera hardware graph Cm,n,l, 
which is an M X grid of l complete bipartite graphs (unit cells). D-Wave 
Two processors use a 512-qubit hardware graph, and the most recent 

processors use a 1152-qubit Ci 2 ,i 2,4 graph. 

The Chimera graph was chosen, in part, because it contains a particularly 
nice clique minor |S| (see Figure [^. This triangle embedding is uniform in the 
sense that each vertex image (or chain) has the same number of vertices, and 
is near-optimal in the sense that it gives a Kj^m minor in Cm,m l, whereas 
Cm,m,l has treewidth LM and therefore contains no Kj^m+2 minoij^ A degree 
argument also shows that any uniform K^m minor requires chains of size at 
least M, while the triangle embedding has chains of size M + 1. 

In practice, a given processor will have a number of inoperable qubits. If 
there are t inoperable qubits, up to t chains in a particular clique embedding 
can be rendered useless. These inoperable qubits force us to find a clique 
minor in an induced subgraph of Cm,n,l- In the face of this, note that there 

^ Taking the triangle embedding and making an image of all the unused qubits gives a 
Klm +1 minor. 
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Fig. 1 (Color online) “Triangle” clique embedding in 04 ^ 4 ^ 4 , which motivated the design of 
the Chimera graph. Every chain in this embedding has 5 qubits. 


are at least four triangle embeddings, which we can find by simply rotating 
the embedding shown in Figure [l] So a first attempt at minimizing the impact 
of inoperable qubits is to choose the triangle embedding for which the greatest 
number of chains survive. 

Triangle embeddings can be generalized further. A triangle embedding con¬ 
sists of overlapping ell-shaped (L-shaped) “bundles” of chains, and each chain 
in a bundle joins a horizontal “wire” with a vertical “wire” via a matching at 
the corner of the bundle. First, the structure of overlapping ell-shaped bundles 
can be generalized from the triangle (we can avoid all four corner unit cells, for 
example). Second, the corner matchings can be chosen arbitrarily to minimize 
the impact of inoperable qubits (if there are three intact vertical wires and 
three intact horizontal wires, we can ensure that they are matched together to 
make three intact chains). These generalizations result in exponential expan¬ 
sion of the number of clique embeddings available, but we can optimize over 
them in polynomial time using a dynamic programming approach. Defining 
and efficiently optimizing over these clique embeddings are the main goals of 
this work. 

In the next section we formalize the definition of native clique embeddings 
that generalize triangle embeddings, and give a combinatorial characterization 
of the same. In Sectionj^we give a dynamic programming technique that, given 
an induced subgraph of Cm,n,l, finds a maximum-sized native clique embed¬ 
ding in polynomial time. One desirable feature of native clique embeddings 
is uniform chain length, which results in uniform, predictable chain dynam¬ 
ics throughout the anneal |13j . Clique embeddings found through heuristic 
methods such as the algorithm described by Cai et al. [3] generally lack this 
property. In Section we compare the results of our approach to those of 
the somewhat similar approach by Klymko et al. m- The shorter chains and 
larger cliques generated by our approach lead to improved tunneling dynamics 
and error insensitivity [Tl lTSllM] . 
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Fig. 2 (Color online) (1) A maximum native clique embedding of K 24 in an induced sub¬ 
graph of C8,8,4 with 26 randomly selected vertices deleted, (r) The corresponding block 
clique embedding, with dots indicating corners of the ell blocks. 


2 Native Clique Embeddings 

We now formally define the structures required to construct and analyze native 
clique embeddings. 

Recall that Chimera Cm,A f.L is an MxN grid of i^L.L unit cells. Specifically, 
Cm,n,l has vertices V = {1, • • • ,M} x {1, • • • , N} x {0,1} x {1, • • • , L} and 
edges: 

{x, y, 0, k) ^ {x + 1, y, 0, k) (horizontal inter-cell couplings), 

{x, y, 1, k) ^ {x,y + 1,1, k) (vertical inter-cell couplings), and 

{x, y, 0, fci) ~ {x, y, 1, k 2 ) (intra-cell couplings). 

We construct native clique embeddings using wires. For t > 1, a horizontal 
wire of length t is a contiguous set of vertices {{x + i,y,0,k) : i G [0, t — 1]}, 
whose induced subgraph is a path on t vertices. Likewise, a vertical wire is a 
set {(x, y -|- i, 1, fc) : z € [0, t — 1]}. An ell is the union of a horizontal wire and 

a vertical wire where there is an edge between one end of the horizontal wire 

and one end of the vertical wire. Note that these ends are necessarily in the 
same unit cell, which we call the corner of the ell; for an ell £ we denote the 
corner by c{i). In the embeddings we study, each chain is an ell. Our aim now 
is to specify an orderly way of arranging them into a clique minor. 

Looking at Figure one may notice that chains appear in sets that inter¬ 
sect the same unit cells. With this in mind, for an ell I we define its ell block 
{X{£),c{£)): X{£) is the set of unit cells intersecting i, and again c{£) is the 
corner of £, which we must specify for the ell block in order to avoid ambiguity 
in the case of horizontal or vertical wires of length 1. We define an ell bundle 
R as a (possibly empty) set of vertex-disjoint ells £i,... ,£p with the same ell 
blocks, i.e. such that |{(X(f), c(f)) | £ G B}\ < 1. 
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A block clique embedding is a set X of n ell blocks {(ATi, ci),..., (Al„, c„)} 
such that 

— each Xi contains n unit cells (so ells have length n + 1), and 

— every distinct pair Xi, Xj in X intersects at exactly one unit cell, which 
is in the horizontal component of one ell block and the vertical component 
of the other. 

A native clique embedding respecting a block clique embedding A is a 
collection B of ell bundles {Bi, ... ,Bn] such that for each i and for each 
£ € Bi, {Xi,Ci) = (X(£), c(£)), i.e. {Xi,Ci) is the ell block for each ell in Bi. 
From this definition and the above, we infer that 

— any two ells £ and £' in the same bundle have exactly two edges between 
them, both in the unit cell c(t') = c(£'), and 

— any two ells £ and £' in different bundles have exactly one edge between 
them, and it is in the unit cell X{£) n X{£''). 

Hence a native clique embedding is a clique embedding. 

It turns out that in a block clique embedding X = {(Ai, ci),..., (A„, c„)}, 
the corners ci,..., form a permutation in the nx n matrix representing the 
unit cells of the graph Cn,n,L- These permutations have a specific structure 
that is in direct correspondence with a class of permutations representable by 
circular point sets studied recently by Vatter and Waton [T^ . 

The following theorem provides a constructive classification of block clique 
embeddings and shows that, in contrast with triangle embeddings, native 
clique embeddings exist in abundance. 

Theorem 1 In a Cn^n.L Chimera graph for n>2, there are 4"“^ block clique 
embeddings that contain n ell blocks. In particular, they are in natural bijeetion 
with the set {E,W} x {NE, NW, SE, SW}”-^ x {N,S}. 

To prove Theorem we first show that each ell block in a block clique 
embedding has a distinct shape. Define the width of an ell as the number of 
vertices in its horizontal wire, and its height as the number of vertices in its 
vertical wire. All ells in an ell bundle have the same width and height, so define 
the width and height of an ell bundle or an ell block as the width and height 
of its constituent ells. 

Lemma 1 Let X = {(Ai, ci),..., (A„, c„)} be a block clique embedding in 
Cn,n,L- Then the ell blocks of X have distinct heights. 

Proof We will establish that X has a unique ell block of height i for each 
1 < i < n. Clearly two ells of height 1 or two ells of height n cannot intersect 
properly (i.e. in exactly one unit which is horizontal for one ell block and 
vertical for the other). Assume for contradiction that X contains two ells (A, c) 
and (A', c') of height 1 < i < n. Their horizontal and vertical components 
must occupy different rows and columns respectively. Up to symmetry, there 
are two cases to consider where ells of the same height intersect properly. 
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Case 1 Case 2 


Fig. 3 Two cases of intersecting ell blocks with the same shape. 


E*S*S* W* 


Fig. 4 Two examples of the first ell block’s orientation and location. * denotes a (possibly 
empty) subword containing zero Ss. 


shown in Figure]^ In both cases, we name the upper ell block (X,c). Since 
there are n ell blocks and n cells per block, every non-corner cell of every ell 
block must intersect another ell block. In Case I, we have shaded cells for which 
the horizontal or vertical coordinate is unique among the two ells. Consider an 
ell block (y, d) which properly intersects (X, c) in a shaded cell. To properly 
intersect (X',c'), (X, c) cannot intersect (X, c) again, so must intersect (X',c') 
in a shaded cell also. Therefore, Y must size greater than n, a contradiction. 
In Case 2, we have shaded the cell which lies directly north of a corner. In this 
case, it is clear that no ell block can properly intersect (X, c) at the gray cell 
and also intersect (X', c') properly. □ 

Without loss of generality, we will assign labels to each ell block (X^, Ci) G 
X so that Xi has height i and width n — i + 1. 

Proof (of Theorem^ We hrst give a mapping from the set of 4”“^ words to 
the set of block clique embeddings and then provide the inverse mapping. 

Let W = si,..., Sn be a word for which si G {E, W}, each of S 2 , • ■ •, Sn-i is 
in {NE, NW, SE, SW}, and Sn G {N, S}. We construct a block clique embedding 
X = {(Xi, Cl),..., (X„, c„)} from W in such a way that each ell block (Xi, a) 
corresponds to the subword s^. We denote the location of a corner Ci by its 
Cartesian coordinates (xi,yi) G {!,...,n}^ and coordinates increase to the 
east (for x) and north (for y). 

If Si = W we place ci so that it is west of the remaining corners, i.e. xi = 1. 
If Si = E, we place ci so that it is east of the remaining corners, i.e. xi = n. 
We select yi so that there are yi — 1 Ss following the subword si in the word 
W. In Figure we show two initial ell block placements for words where * 
denotes a (possibly empty) subword containing zero Ss. 

The placement of the first ell block defines a working rectangle i?i of cells 
intersecting horizontal wires but not vertical wires in the partially-constructed 
block clique embedding. Every time we place an ell block, the working rect¬ 
angle becomes one unit taller and narrower, and the vertical component of 
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£••• ENE--- ENENW--- ENENWSE--- ENENWSES 


Fig. 5 An example of the iterative construction of a block clique embedding. 


the recently-placed ell covers either the left-most or right-most column of the 
previous working rectangle. 

Suppose the working rectangle after the selection of i ell blocks is Ri = 
[x,x'] X [y,y'] in Cartesian coordinates, where x < x' and y < y'. We choose 
corner Ci based on subword Si as follows, noting that x' = x + n — i — \ and 

y' =y + i-l: 

— If Si = NE or N, a = {x',y' + 1). 

— If Si = NW or N, Ci = {x, y' + 1). 

— If Si = SE or S, Ci = (a;', y — 1). 

— If Si = SW or S, Ci = {x,y — I). 

Letting a = {xi,yi), we then choose Xi = ({cci} x [y,y']) U {[x,x’] x {yi}) and 
update the working rectangle i?i+i = ([x,a:']\{a;i}) x ([?/,?/'] U {yi}). Observe 
that every time that we find a W, x increases by 1 and every time we find a E, 
x' decreases by 1; every time we find an N, y increases by 1 and every time an 
S, y' decrease by 1. Figure shows an example of how the construction might 
proceed. 

Note that Ri never shares a column with a corner in {cj | j < *}, and 
the set of rows intersected by Ri is the same set of rows intersected by 
{cj I j < *}. It therefore follows from the construction that {ci,...,c„} is 
a permutation. Furthermore, by the construction, for i < j there is always a 
unit cell Xi n Xj in the working rectangle Ri - specifically, it is {xj,yi). So 
X = {(Vi, Cl),..., (V„, c„)} is indeed a block clique embedding. 

Now, we invert our construction to show that it is in fact bijective. Let 
X = {(Xi, Cl),..., (V„, c„)} be a block clique embedding in which Xi has 
height i; we will reconstruct our word W. We call {xn,yi) = Xi n X„ the 
center of X, and say that a point (a:, y) lies to the east or west of the center if 
X > Xn or X < Xn respectively, and likewise north and south. 

— Let Si = E if the corner ci lies to the east of the center, and W otherwise. 

— For 1 < f < n, let Si = Oifoi, where Oi = N if Ci lies to the north of the 
center, and S otherwise; and 6i = E if Ci lies to the east of the center, and 
W otherwise. 

— Let s„ = N if the corner c„ lies to the north of the center, and S otherwise. 

Therefore, W = S 1 S 2 ... s„ is in {E, W} x {NE, NW, SE, SW}"“^ x (N, S}. To 
see that this construction is in fact the inverse, we observe that each ell block 
(Xi,Ci) is associated to the direction Si in which the corner Ci lies from the 
center, so we obtain the same word we began with. □ 
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Fig. 6 The four rectangles Ri = Rto{Xi,Ci) with (Xi,Ci) S Xto{R)- Gray ell blocks repre¬ 
sent precomputed maxPartialEmbedding(Ri) as in Lemmaj^ 


3 Finding optimal native clique embeddings in induced subgraphs 

In this section we describe an algorithm to find a largest native clique embed¬ 
ding B = {Bi ^..., Bn\ in an induced subgraph G of Cm,n,l with n < M, N. 
Our algorithm necessarily takes as input a parameter n, which determines the 
size of the chains, i.e. n + 1. 

Our algorithm is informed by the proof of Theorem We use dynamic 
programming to maximize the block clique embeddings with each working 
rectangle i?, and do so in an orderly way which results in a polynomial-time 
algorithm. 

As a preprocessing step, for each ell block {X, c) we compute a maximum 
bundle and store it as maxBundle{X, c) (it is straightforward to do this in 
0{nL) time per ell block). We then use this information to construct block 
clique embeddings as in the proof of Theorem adding one ell block at a 
time, with working rectangles of increasing height. Let (X,c) be an ell block 
of height L If 1 < z < n — 1, there is a unique working rectangle Rfrom{X,c) 
that can be in effect immediately after {X, c) is placed. If 2 < z < n, there is 
a unique working rectangle Rto [X, c) that can be in effect immediately before 
{X, c) is placed. Note that for each working rectangle i?, the sets 

Xfrom{R) ■■= {{X,C)\R= Rto{X,c)} 


and 

XUR) := {{X,c)\R = Rfroin{X,c)} 

each have size at most four (see Figure]^. 

For any a set of ell bundles B = {Bi, • • • , B„} where each Bi is contained 
in the ell block (W, a) and X = {(Aii, ci), • • • , (X„, c„)}. 




n 

maxBundle{Xi,Ci) 

i=l 


> 


i=l 


provided that the ell blocks are all distinct. This enables us to construct max¬ 
imum native clique embeddings while only considering the ell blocks involved. 

For each working rectangle R of height z (and width n — i), our algorithm 
computes and stores a partial block clique embedding with that working rect¬ 
angle: a set of ell blocks A, = {(ATi, ci), • • • , (W) g)} such that 


R — Rfrom {Xi , Ci ) 
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— Xj has height j for all 1 < j i 

, Cj_ 1 ) for all 2 ^ J ^ h 

In particular, we compute maximum partial block clique embeddings with 
working rectangle i?; those which maximize ||Ai|| over all partial block clique 
embeddings with a given working rectangle. We denote a particular max¬ 
imum partial block clique embedding for a given working rectangle R by 
maxPartialEmbedding{R), though many may exist. Our algorithm will operate 
by extending maximum partial block clique embeddings by ells, so we define 
maxExtension[X, c) = maxPartialEmbedding{Rto{X, c)) U {(X, c)} when X 
has height h> 1, and maxExtension{X, c) = {(X, c)} otherwise. 

The following lemma encodes the key step of the algorithm: to find a max¬ 
imum partial block clique embedding with working rectangle R, we need only 
consider partial block clique embeddings that include an ell block (X, c) such 
that R = Rfrom{X,c). 

Lemma 2 Given a rectangle R of height i > 1 and width n — i, 

\\maxPartialEmbedding{R)\\ = max \maxExtension(X, c)\. 

(x,c)eXt„{R) 

Proof We proceed by induction on i. When i = 1, our claim follows imme¬ 
diately from definitions. Assume that i > 1 and that our claim holds for all 
working rectangles of height i — 1. 

We consider the complete set of partial block clique embeddings with 
working rectangle R, S = {Xi \ Rfrom.{Xi,c,) = R} = {Xi \ (Xi,c,) G 
Xto{R^)} where {Xi,Ci) is the ell block with height i in Xi. By definition, 
\\maxPartialEmbedding{R)\\ = max;t’ieS IIAi||. For contradiction, pick some 
maximum Xi G S and suppose that 

||Ai|| > \\maxPartialEmhedding{Rto{X,c))\J {{X,c)'\\\ 
for all (X, c) £ Xto{R). In particular, letting Xi_i = Xi \ {{Xi, Ci)}, 

||X,|| = ||X,_i|l + ||{(X„c.)}|| 

> \\maxPartialEmbedding{RtoiXi,Ci))\\ + \maxBundle{Xi,Ci)\, 
a contradiction since Xi-i has working rectangle Rto{Xi, Ci). □ 

We now present the algorithm. The idea follows from Lemma to com¬ 
pute maxPartialEmbedding{R) for rectangles of increasing height. To do so 
we treat the set of possible working rectangles as a digraph, where i? —>■ i?' 
if and only if there is an ell block (X,c) for which R G Rfrom{X,c) and 
R' G Rto{X,c). This means that if i? — R', the height of R' is one more 
than the height of R. The number of edges in this digraph is equal to the 
number of ell blocks. To compute maxPartialEmbedding{R'), assuming that 
we have computed maxPartialEmbedding{R) for all rectangles of lesser height, 
we simply set maxPartialEmbedding{R') to be a maximum partial block clique 
embedding in the set 

{maxExtension{X, c) \ (X, c) G Xto(R')} ■ 
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Once we have computed maxPartialEmbedding{R) for all rectangles R of height 
n — 1 and width 1, we pick a maximum-sized clique embedding from the set 

{maxExtension{X, c) \ {X, c) has height n — 1 and width 1} . 

Pseudocode is given in Algorithm 


Algorithm 1 The algorithm to find a maximum-sized native clique embedding 
in an induced subgraph of a Chimera graph. 

1: function NativeCliqueEmbed(G, n) 

2: for i = 1,..., n — 1 do 

3: for each rectangle R of height i and width n — i do 

4: maxPartialEmbedding{R) <— 0 

5: for each ell block {X, c) of height i and width n — i + 1 do 

6: B <— maxPartialEmbedding{Rto{X,c)) U {(X, c)} 

7: if \\maxPartialEmbedding{Rf^g^{X,c))\\ < ||B|| then 

8: maxPartialEmbedding{Rfr„^{X, c)) <— B 

9: Bmax ■<— 0 

10: for each ell block {X, c) of height n and width 1 do 

11: B <— maxPartialEmbedding(Rto{X,c)) U {(X, c)} 

12: if IIBmaxll < ||B|| then 

13: Bmax ^ B 

14: return {maxBundle{X,c,G) \ {X,c) £ Bmax} 


Theorem 2 The NativeCliqueEmbed algorithm finds a maximum-sized native 
clique embedding with chain length n -\-1 in polynomial time. 


Proof We first prove correctness, then the bound on running time. 

Note that the loop beginning on line of Algorithm iterates over all ell 
blocks of height i and width n — i. Given a rectangle i?, there are up to four ell 
blocks (X, c) for which R = Xto (X, c). If we ignore all ell blocks except those 
incident to a particular rectangle R, this loop implements Lemma directly. 
Therefore, when we reach line|^ we have computed maxPartialEmbedding{R) 
for all R with height n — 1 and width 1. 

Let B = {Bi, ■ ■ ■ ,Bn} be a maxmimum native clique embedding where 
Bi is an ell bundle in the ell block (Xi,Ci) with height i. By Theorem]^ 
Xn-i = {(Xi, Cl), • • • , (X„_i, Cn-i)} is a partial block clique embedding with 
working rectangle RtoiXn, c„). By Lemma 

\\maxPartialEmbedding{RtoiXn, Cn))\\ > HAn-iH , 


so we see a clique embedding of size at least HAn-iH -1- \maxBundle{Xn,Cn)\ 
when line 11 is reached with (X, c) = (X„,c„). Therefore, when line [Td] is 
reached, USmosll > ll^ll and a native clique embedding of maximum size has 
been found. 

We can compute maxBundle{X, c, G) in 0{nL) time, and there are poly- 
nomially many ell blocks and rectangles: There are at most MN possible lo¬ 
cations of a rectangle’s lower-left corner, and n possible shapes, thus at most 
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nMN rectangles. Likewise there are at most MN possible locations for an ell 
block’s corner, and at most 4(n — 1) ell blocks containing n unit cells with a 
given corner, thus at most AnMN ell blocks. 

It follows that each line in Algorithm is evaluated 0{nMN) times, and 
the preprocessing step of computing maxBundle{X,c^G) for each ell block 
{X,c) naively takes 0{'n?MNL) time. Consequently, with the rough bound 
that each line in Algorithm takes OiriL) time for a single evaluation, we can 
bound the total running time of our algorithm by 0{n?MNL). □ 

Remark. The 0{n^MNL)-tim.e bound on Algorithm[^is quadratic in the num¬ 
ber of vertices in Cn,m,l, i.e. 0{'n?MNL) C 0{{MNLY). Assume M < N 
and L is constant. With a little more care, we can modify Algorithm to 
achieve a bound of 0{N^) instead of 0{N'^). Doing this involves (a) precom¬ 
puting all maximum horizontal and vertical line bundles in time 0{LNM‘^ + 
LIVIN'^) with a dynamic programming approach, which allows us to compute 
\maxBundle{X, c, G)! in 0(1) time, and (b) exploiting the fact that throughout 
the algorithm, we need only keep track of the size of maximum partial embed¬ 
dings and the route used to reach it (replacing maxPartialEmbedding{R) with 
a mapping R i—(V, c) S Xto{R)), rather than the embeddings themselves. 

NativeCliqueEmbed gives a maximum native clique embedding for a fixed 
chain length. To find a maximum native clique embedding over all chain 
lengths for a given graph, we simply repeat the process for each choice of 
n G {2,... ,min{M, V}}. For M < N and constant L, this gives an overall 
running time on a subgraph of Cm,n,l of 0{N^) with the naive implementation 
and O(V^) with the refinement discussed above. 


3.1 Induced and general subgraphs 

We now discuss the motivation of using induced subgraphs and how to ap¬ 
proach more general subgraphs. 

Recall that we have restricted our attention to induced subgraphs rather 
than more general subgraphs because failed couplers adjoining working qubits 
are relatively rare. In an induced subgraph, we will still focus on horizontal 
and vertical wires, and it is easy to find a maximum set of such wires in a line 
of unit cells. 

A maximum ell bundle is a maximum-sized set of vertex-disjoint ells that 
occupy an ell block, and the size of an ell bundle is the number of ells contained 
therein. It is simple to find a maximum ell bundle in a given ell block: finding 
maximum sets Sr and Sy of horizontal and vertical wires spanning a line of 
cells is trivial, and these lines may be paired off arbitrarily since the corner 
is a complete bipartite graph. Here there may be difficulty in generalizing 
even this relatively simple optimization problem in the face of arbitrary edge 
deletion. Finding a maximum ell bundle in an arbitrary Chimera subgraph is 
polynomially equivalent to finding a maximum clique in D 2 {C{B)) where B 
is a bipartite graph, C{B) is the line graph of B, and D 2 {G) is the distance-2 
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graph with vertex set V{G) and edge set {uv : da{u,v) = 2}. We are unsure 
of the complexity of this problem, but expect that it is NP-complete. 

Note that we can easily relax the requirement of an induced subgraph when 
couplers between unit cells are defective - in any case, one just computes the 
number of wires in a line of cells where all qubits and couplers are contained 
in the subgraph. In short, failed inter-cell couplers don’t increase the difficulty 
of the problem. 

If we restrict our attention to Chimera graphs Cm,n,l where L = 0(log TV), 
or assume L to be a constant, then this difficulty at the corners can be swept 
under the rug. At present, this appears to be a reasonable consideration, as 
it is much easier from a manufacturing and design standpoint to increase N 
than it is to increase L. 

However, there is a further challenge introduced by intra-cell couplers. Our 
algorithm intrinsically relies upon the assumption that couplers exist between 
any two ells whose unit cells intersect. Missing intra-cell couplers void that 
assumption. The easiest remedy for this obstruction is to consider vertex covers 
of the failed intra-cell couplers: Given a graph G that is the subgraph of 
Cm,n,l induced by vertex set W, i.e. G = Cm,n,l[W], let Ui,... ,Us be the 
list of minimal vertex covers of the failed intra-cell couplers. Then, for each 
1 < i < s, compute a largest clique minor considering G\W] for the purpose 
of constructing ells, and G\W \ Ui] for the purpose of growing cliques. This 
approach gives a fixed-parameter tractable algorithm for finding the maximum 
native clique embedding in an arbitrary subgraph of Chimera in terms of the 
number of missing edges with both endpoints intact - for the D-Wave Two 
processors installed at NASA Ames [13] and ISI [T|, this parameter was zero. 


4 Comparison with previous work 

Klymko, Sullivan and Humble gave a greedy embedding algorithm that quickly 
produces similar embeddings to those in this paper nnj. Their algorithm pro¬ 
duces plus-shaped chains (having nearly twice as many qubits), and restricts its 
search to embeddings where the vertical and horizontal components of pluses 
are only allowed to meet at the diagonal of a fixed square. While our algo¬ 
rithm is slower, taking 0{N'^) time compared to their 0{N^) for Cn,n, 4 , it is 
exhaustive, empirically embeds larger cliques, and produces chains of roughly 
half the size. 

Given a family of clique minors, the clique yield of a graph G over that 
family is the size of the largest clique minor in the family, in G. In Figure 
we compare ell- and plus-shaped chains as the grid size grows for several fixed 
percentages of operational qubits. In both families, a similar asymptotic be¬ 
havior becomes clear: for a fixed qubit failure rate, increasing the grid size 
gives diminishing returns in terms of clique yield. However, the difference be¬ 
tween these curves is significant, with ell-shaped chains producing much larger 
clique minors. 



Fast clique minor generation in Chimera qubit connectivity graphs 
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Fig. 7 (Color online) Comparison of clique yields for ell-shaped (red) and plus-shaped 
(blue) chains. The solid lines denote the median, and the shaded regions encompass the 
middle two quartiles. 
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